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INTRODUCTION 

In some engineering applications such as chaotic encryption, chaotic maps have to 
exhibit required statistical and spectral properties close to those of random signals. There 
is a gro wing industrial interest to consider and study thoroughly the property of such map 

Very often, dynamical systems in several dimensions are obtained coupling 1- 
dimension ones and their properties are strongly linked Q5D . 

Quasi-periodic or chaotic motion is frequently present in complicated dynamical 
systems whereas simple dynamical systems often involve only periodic motion. The 
most famous theorem in this field of research is the Sharkovskii's theorem, which 
addresses the existence of periodic orbits of continuous maps of the real line into itself. 
This theorem was once proved toward the year 1962 and published only two years after 

ill- 
Mathematical results concerning periodic orbits are often obtained for functions on 
real intervals. However, most of the time, as the complex behaviour of chaotic dynami- 
cal systems is not explicitly tractable, mathematicians have recourse to computer simu- 
lations. The main question which arises then is: does these numerical computations are 
reliable ? 

As an example we report the results of some computer experiments on the orbit 
structure of the discrete maps on a finite set which arise when the logistic map is iterated 
"naively" on the computer. 



Due to the discrete nature of floating points used by computers, there is a huge gap 
between these results and the theoretical results obtained when this map is considered 
on a real interval. This gap can be narrowed in some sense (i.e. avoidingthe collapse of 
periodic orbits) in higher dimensions when ultra weak coupling is used flsLZD. 

Nowadays the claim is to understand precisely which periodic orbit can be observed 
numerically in such systems. In a first attempt we study in this paper the orbits generated 
by the iterations of a one-dimensional system on a finite set with a cardinal N. The 
final goal of a good understanding of the actual behaviour of dynamical systems acting 
on floating numbers (i.e. the numbers used by computers) will be only reached after this 
first step will be achieved. 

On finite set, only periodic orbits can exist. For a given function we can compute all 
the orbits, all together they form a global orbit pattern. We formalise such a gop as the 
ordered set of periods when the initial value thumbs the finite set in the increasing order. 
We are able to predict, using closed formulas, the number of gop for the set J^)v of all 
the functions on X. We also explore by computer experiments special subsets of J^jy> 
such as sets of locally "rigid" functions which presents interesting patterns of gop. 

This article is organized as follows : in the section "Computational divergences" we 
display some examples of such computational divergences for the logistic map in various 
ways of discretization. In the section "Pattern defined by all the orbits of a dynamical 
system" we introduce a new mathematical tool: the global orbit pattern, in order to 
describe more precisely the behaviour of dynamical systems on finite sets. In the section 
"Cardinal of classes" we give some closed formulas related to the cardinal of classes of 
gop of J£}y. In the section "Functions with local properties" we study the case of sets of 
functions with a kind of local "rigidity" versus their gop, in order to show the usefulness 
of these new tools. 



COMPUTATIONAL DIVERGENCES 

Discretized logistic map 

As an example of collapsing effects which happen when using computers in numerical 
experiments, we presents the results of a sampling study in double precision of a 
discretization of the logistic map fa : [0, 1] — > [0, 1] (see Fig.© 

f 4 (x)=4x{l-x) (1) 
and its associated dynamical system 

x n+ i=4x n (l-x n ) (2) 

which has excellent ergodic properties on the real interval. 
There exists an unstable fixed point 0. 

The set j 1 ^, ^fi} = {0.3454915,0.9045084} is the period-2 orbit. 

In fact there exist an infinity of periodic orbits and an infinity of periods. This dynamical 
system possesses an invariant measure (see Fig. [2]): 



FIGURE 1. Graph of the map f(x) = 4x(l -x) on [0, 1] 
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FIGURE 2. Invariant measure of the logistic map 



P{*) = m c (3) 

lZyfX( \ — x) 

However, in numerical computations using ordinary (IEEE-754) double precision 
numbers - so that the working interval contains of the order of 10 16 representable points 
- out of 1,000 randomly chosen initial points (see Tabled)), 

• 596, i.e., the majority, converged to the fixed point corresponding to the unstable 
fixed point {0} in equation |2l 

• 404 converged to a cycle of period 15,784,521. 

Thus, in this case at least, the very long-term behaviour of numerical orbits is, for a 
substantial fraction of initial points, in flagrant disagreement with the true behaviour of 
typical orbits of the original smooth logistic map. 

In others numerical experiments we have performed, the computer working with fixed 
finite precision is able to represent finitely many points in the interval in question. It is 



probably good, for purposes of orientation, to think of the case where the representable 
points are uniformly spaced in the interval. The true logistic map is then approximated 
by a discretized map, sending the finite set of representable points in the interval to itself. 

Describing the discretized mapping exactly is usually complicated, but it is roughly 
the mapping obtained by applying the exact smooth mapping to each of the discrete 
representable points and "rounding" the result to the nearest representable point. In our 
experiments uniformly spaced points in the interval with several order of discretiza- 
tion (ranging from 9 to 2,001 points) are involved. In each experiment the questions 
addressed are: 

• how many periodic cycles are there and what are their periods ? 

• how large are their respective basins of attraction, i.e. , for each periodic cycle, how 
many initial points give orbits with eventually land on the cycle in question ? 



TABLE 1. Coexisting periodic orbits found using 1 ,000 random 
initial points for double precision numbers 



Period 


Orbit 


Relative Basin size 


1 

15,784,521 


{0} (unstable fixed point) 
Scattered over the interval 


596 over 1,000 
404 over 1,000 



TABLE 2. Coexisting periodic orbits 
for the discretization with regular meshes 



of JV 


= 9, 10 and 1 1 points 




N 


Period 


Orbit 


Basin size 


9 


1 


{0} 


3 over 9 


9 


1 


{6} 


2 over 9 


9 


1 


{3,7} 


4 over 9 


10 


1 


{0} 


2 over 10 


10 


2 


{3,8} 


8 over 10 


11 


1 


{0} 


3 over 1 1 


11 


4 


{3,8,6,9} 


8 over 1 1 



TABLE 3. Coexisting periodic orbits for the discretization with 
regular meshes of N = 99, 100 and 101 points 



N 


Period 


Orbit 


Basin size 


99 


1 


{0} 


3 over 99 


99 


10 


{3,11,39,93,18,58,94,15,50,97} 


96 over 99 


100 


1 


{0} 


2 over 100 


100 


1 


{74} 


2 over 100 


100 


6 


{11,39,94,18,58,96} 


72 over 100 


100 


7 


{7,26,76,70,82,56,97} 


24 over 100 


101 


1 


{0} 


3 over 101 


101 


1 


{75} 


2 over 101 


101 


1 


{16,61,95} 


96 over 101 



On an another hand, for relatively coarse discretizations the orbit structure is deter- 
mined completely, i.e., all the periodic cycles and the exact sizes of their basins of attrac- 
tion are found. Some representative results are given in Tables [2] to HI In theses tables, N 



TABLE 4. Coexisting periodic orbits for the discretization with regular meshes 
of jV = 1,999;2,000 and 2,001 points 



N 


Period 


Orbit 


Basin size 


1,999 


i 


{0} 


3 over 1,999 


1,999 


4 


{554; 1,601; 1,272; 1,848} 


990 over 1,999 


1,999 


8 


{3; 11;43; 168; 615; 1,702; 1,008; 1,997} 


1,006 over 1,999 


2,000 


1 


{0} 


2 over 2,000 


2,000 


1 


{1,499} 


1 a ^ r\r\r\ 

14 over 2,000 


2,000 


2 


{691; 1,808} 


138 over 2,000 


2,000 


3 


{276; 1,221; 1,900} 


6 over 2,000 


2,000 


8 


{3; 11;43;168;615; 1,703; 1,008; 1,998} 


1,840 over 2,000 


2,001 


1 


{0} 


5 over 2,001 


2,001 


1 


{1,500} 


34 over 2,001 


2,001 


2 


{691; 1,809} 


92 over 2,001 


2,001 


8 


{3; 11;43;168;615; 1,703; 1,011; 1,999} 


608 over 2,001 


2,001 


18 


{35; 137;510; 1,519; 1,461; 1,574;...} 


263 over 2,001 


2,001 


25 


{27;106;401;1,282;1,840;588;...} 


1,262 over 2,001 



denotes the order of the discretization, i.e., the representable points are the numbers, -p, 
with < j < N. 

The Table[2]shows coexisting periodic orbits for the discretization with regular meshes 
of N = 9, 10 and 1 1 points. There are exactly 3, 2 and 2 cycles. 

The Table [3] shows coexisting periodic orbits for the discretization with regular meshes 
ofN = 99, 100 and 101 points. There are exactly 2, 4 and 3 cycles. 
The Table |4] shows coexisting periodic orbits for the discretization with regular meshes 
ofN= 1, 999,^ = 2,000 and N = 2, 001 points. 

It seems that the computation of numerical approximations of the periodic orbits leads 
to unpredictable results. 



Statistical properties 

Many others examples could be given, but those given may serve to illustrate the 
intriguing character of the results: the outcomes proves to be extremely sensitive to the 
details of the experiment, but the results all have a similar flavour : a relatively small 
number of cycles attract near all orbits, and the lengths of these significant cycles are 
much larger than one but much smaller than the number of representable points. 

In 01, P. Diamond and A. Pokrovskii, suggest that statistical properties of the phe- 
nomenon of computational collapse of discretized chaotic mapping can be modelled by 
random mappings with an absorbing centre. The model gives results which are very 
much in line with computational experiments and there appears to be a type of univer- 
sality summarised by an Arcsine law. The effects are discussed with special reference to 
the family of mappings 

x n+l = 1 - |1 - 2x n \ £ < x < 1 1 < i < 2 (4) 
Computer experiments show close agreement with prediction of the model. 



However these results are of statistical nature, they do not give accurate information 
on the exact nature of the orbits (e.g. length of the shortest one, of the greater one, 
size of their basin of attraction ...). It is why we consider the problem of computational 
discrepancies in an original way in the next section. 

PATTERN DEFINED BY ALL THE ORBITS OF A DYNAMICAL 

SYSTEM 

In this section in order to describe more precisely which kind of behaviour occurs in 
discretized dynamical systems on finite sets we conceive a new mathematical tool: the 
global orbit pattern of a function that is the set of the periods of every different orbits 
of the dynamical system associated to the function when the initial points are took in 
increasing order. 

General definitions 

For every xq EX, let {x,} be the sequence of the orbit of the dynamical system 
associated to the function / which maps X onto X defined by 

Xi+i = f(xi) for i > 0. (5) 
For convenience \/xo G X we denote 

f(x )=x (6) 

and 

Vp > l,Vx G X, fP(x ) = /o/o...o/ (x ). (7) 

P times 

Hence 

Xi = f(x ). (8) 

The orbit of xq under/is the set of points &(xo,f) = {f l (xo),i > 0} = {xi,i > 0}. 
The starting point xq for the orbit is called the initial value of the orbit. 
A point x is a fixed point of the map / if f(x) = x. 

A point x is a periodic point with period p if f p (x) = x and f k (x) ^ x for all k such that 
< k < p, p is called the order of x. 

If x is periodic of order p, then the orbit of x under / is the finite set 

{x,f(x),f 2 {x), . . . ,f p ~ l (x)}. We will call this set the periodic orbit of order p or a 

/>-cycle. 

A fixed point is then a 1 -cycle. 

The point x is an eventually periodic point of / with order p if there exists K > such 
that Vk > K f k+ P (x) = f (x) . 

G X, we denote G){x,f) the order of x under / or simply co(x) when the map / 
involved is obvious. 



A subset T of X is invariant under / if / (T) = T. That is equivalent to say that T is 
invariant under / if and only if f(T) C T and f~ l (T) C T. 

Notation pC is the cardinal of the finite set X. 



Map on finite set 

Along this paper, N is a non-zero integer and (JA stands for the cardinal of any 
finite set A. In this article, we consider X as an ordered finite set with N elements. 
We denote it X N , it is isomorphic to the interval [[0,iV- 1]] C N. Then fJXjy = N. Let 
/ be a map from X^ into X^. We denote by J^n the set of the maps from X^ into 
Xn. Clearly, is a finite set and ftJ^v = N N elements. For all x £ Xjv, 0(x,f) is 
necessarily a finite set with at most Af elements. Indeed, let us consider the sequence 
{x,f(x),f 2 (x),...,f N ~ 1 (x),f N (x)} of the first N+\ iterated points. Thanks to the 
Dirichlet's box principle, two elements are equals because X^ has exactly N different 
values. Thus, every initial value of X^ leads ultimately to a repeating cycle. More 
precisely, if x is a fixed point 0(x,f) is the unique element x and if x is a periodic 
point with order p, &(x,f) has exactly p elements. In this case, the orbit of x under / is 
the set @{x,f) = {x 1 f(x) 1 f 2 (x) 1 . . .,f p ^ 1 (x)}. If x is an eventually periodic point with 
order p, there exists K > such that Vk>K f k+p (x) = f (x) . In this case, the orbit of x 
under / is the set 0{xj) = {x,f(x),f 2 (x),. . . ,f K (x),f K+l (x), . . .J k +p-\x)}. 



Equivalence classes 

Components 

Let / £ We consider on X^ the relation ~ defined by : Vx,y £ Xjv, x ~ 
3^ £ N such that / (x) £ &(x' ,f). The relation ~ is an equivalence relation on Xff. 
S^n/ ~ is the collection of the equivalence classes that we will call components of Xpj 
under / which constitute a partition of X^. The number of components are given in 

Asymptotic properties of the number of cycles and components are studied in 
For each component, we take as class representative element the least element of the 
component. The components will be written T^(xq, /),... , T^(x PfN ,f) where Xi is the 
least element of T^(xi : f). 

By analogy with real dynamical systems, we can define attractive and repulsive 
components in discretized dynamical systems as follows. 

Definition 1 A component is repulsive when it is a cycle. Otherwise, the component 
is attractive. 



Remark In other words, a component is attractive when the component contains at 
least an eventually periodic element. The corresponding cycle is strictly contained in an 
attractive component. 

Examples are given in Tables [5l [6] and |7J 

For instance, in Tabled the fonction / has {2,7} as period-2 orbit and {1,2,7,9} as 
component which is attractive because 1 and 2 are eventually periodic elements. 



TABLE 5. Orbits and components of a function belonging to &\ i with gop 



[2,2,1,3]!!. 


Function 




orbit/component/nature 




-> 


6 






1 -> 


3 






2 -> 


2 


period-2 orbit : {6,9} {0,6,9} 


attractive 


3 -> 


5 






4 -> 


8 


period-2 orbit: {5,10} {1,3,5,10} 


attractive 


5 -> 


10 






6 -> 


9 


fixed point : {2} {2} 


repulsive 


7 


4 






8 


7 


period-3 orbit : {4,8,7} {4,8,7} 


repulsive 


9 


6 






10 -> 


5 






TABLE 6. 


Orbits and components of a function belonging to &\ \ with gop 


[2,2,l,3]n. 








Function 




orbit/component/nature 




-> 


4 






1 -> 


2 






2 -> 


7 


period-2 orbit: {4,8} {0,4,8} 


attractive 


3 -> 


3 






4 -> 


8 


period-2 orbit: {2,7} {1,2,7,9} 


attractive 


5 -> 


10 






6 -> 


5 


fixed point : {3} {3} 


repulsive 


7 -> 


2 






8 -> 


4 


period-3 orbit : {5, 10,6} {5, 10,6} 


repulsive 


9 -> 


1 






10 -> 


6 







Order of elements 
Here are some remarks on the order of elements of components. 
Remark The order of every element of a component is the length of its inner cycle. 



TABLE 7. Orbits and components of a function belonging to & 



\ i with gop 



[2,2,1, 3] n . 



Function orbit/component/nature 



0^9 
1 -> 6 



2 -> 4 period-2 orbit: {0,9} {0,9} n 
3^7 

4 -> 10 period-2 orbit: {1,6} {1,6} n 
5^3 

6 -> 1 fixed point: {10} {2,4,8,10} a 
7^5 

8^2 period-3 orbit: {3,7,5} {3,7,5} n 



repulsive 



repulsive 



repulsive 



attractive 



9^0 
10 10 



Definition 2 For all x E X^, there exists i E [[0,p/yv]] such that x belongs to the 
component Tn (•*;,/)• Then co(x,f) is equal to the order (a(xj,f). 

Remark For all i E [[0,p/yv]], TV (*,•,/) is an invariant subset of X N under /. 

In the example given in Table [5] the order of the element is 2, the order of the 
element 1 is 2, the order of the element 4 is 3. The elements 1 and 3 have the same order. 



For each / E J^v, we can determine the components of Xpj under /. For each compo- 
nent, we determine the order of any element. Thus, for each / E J^jv, we have a set of 
orders that we will denote Q,(f,N). Be given /, there exist pt^ components and pf^ 
representative elements such that xq < x\ < . . . < x PfN . 

For each / E ^v, the sequence [(d(xo),(o(xi), . . . ,(o(x PfN );f]j? N with xq < x\ < 
... < x PfN will design the global orbit pattern of / E ^y. 

We will write gop(f) = [co(*o), fi>(*i), . . . , (o(x Pf N );f}^ N . 

When the reference to f E is obvious, we will write shortly 

gop(f) = [G)(x ),(D(xi), ... ,(D(x PfN )] N or gop(f) = [(o(xo),(o(x 1 ),...,(o(x PfN )} . 

For example, the same gop associated to the functions given in Tables [5] [6] and [7] is 



[2,2, 1,3] ii. 

Another example is given in Tabled] In that example, we have 0)(0) = 2, ft>(3) = 1, 
co(4) = 4. 



Definition of global orbit pattern 



TABLE 8. Orbits and components of a function belonging to J?g with 

gop[2,l,4] g . 



Function orbit/component/nature 



- 


-> 1 








1 - 


-> 


period-2 orbit : {0, 1} 


{0,1,2} 


attractive 


2 - 


-> 








3 - 


-> 3 


fixed point : {3} 


{3} 


repulsive 


4 - 


-> 5 








5 - 


-> 6 


period-4 orbit : {4,5,6,7} 


{4,5,6,7} 


repulsive 


6 - 


4 7 








7 - 


■+ 4 









Definition 3 The set of all the global orbit patterns of J^jv is called ( JF#) . 



For example, for N = 5, the set ^(#5) is 
{[1];[1,1];[1,1,1];[1,2];[1,1,1,1];[1,1,2];[1,3];[1,1,1,1,1];[1,1,1,2];[1, 1,2,1]; 
;[1,2,1,1];[1,2,2];[1,3,1];[1,4];;[2,1];[2,1,1];[2,2];[2,1,1,1];[2,1,2];[2,2,1] 
;[3];[3,1];[3,1,1];[3,2];[4];[4,1];[5]}. 



Class of gop 

We give the following definitions : 

Definition 4 Let be A = [ft>i, . . . , G) p ]ff a gop. Then the class of A, written A, is the 
set of all the functions / G &m such that the global orbit pattern associated to / is A. 

For example, for = 11, the class of the gop [2, 2, 1, 3] n contains the following few 
of many functions defined in Tables[5l[6]and|7J The periodic orbit which are encountered 
have the same length nevertheless there are different. 

Definition 5 Let be A = [ft>i , . . . , a} p ]^ a gop. 

p 

Then the modulus of A is \A\ = £ (Ok- 

k=l 

Remark \[(Oi , . . . , co p ]n\ < N. 

Notation [ft^jjv means [ co, ... , coJ n and [ft£, v~\n means [ co, . . . , ft) , v, . .. , vJ n- 

k times k times m times 



Threshold functions 

Ordering the discrete maps 
Theorem 1 The sets J% and |1,A^]] are isomorphic. 

Proof We define the function from J^X? to [[1,A^^]] by : for each / e J^iv, is 
the integer n such that n= £ f(k)N N ^ l ^ k + 1. 
Then is well defined because n e [[1,^]]. 

Let n be a given integer between 1 and A^. We convert n — 1 in base Af : there 
exists a unique Af-tuple (a„_i ; o;«n-i,i; • • • ;«n-i,jv-i) £ [[0,Af — 1]]^ such that 

« — 1 = a n-\ N-l-iN N ~'~ . We can thus define the map f n with : Vz E X^, 

«=o 

/„(/) = a n -\fl-i-\. Then is one to one. 

Remark This implies ^ is totally ordered. 
Definition 6 Let / e =^v. Then 

n = \f{k)N N - x - k + \ (9) 

is called the rank of /. 



Threshold functions 

Be given a global orbit pattern A, we are exploring the class A. 

Theorem 2 For every A e £f(^iv), the class A has a unique function with minimal 
rank. 

Definition 7 For every class A e £f(J^jv), the function defined by the previous 
theorem will be called the threshold function for the class A and will be denoted by 
TV (A) or 7Y(A). 

To prove the theorem, we need the following definition : 

Definition 8 Let / e ^/v be a function. Let p a non zero integer smaller than N. 
Let be xi , . . . , x p p consecutive elements of Xpj. Then x\ , . . . ,x p is a canonical p-cycle in 
relation to / if Vj e - 1]] , f(xj) = x j+ \ and f(x p ) = x\ . 



Proof Let [©i, . . . , a> p ] be a global orbit pattern of We construct a specific 

function / belonging to the class [a>i , . . . , © p ] and we prove that the function so obtained 
is the smallest with respect to the order on J^}y. With the first ©i elements of [[0,iV— 1]], 
that is the set of integers [[0, ©i — 1]], we construct the canonical ©i -cycle : if C0\ = 1, 
we define f(0) = 0, else /(0) = 1, /(l) = 2, . .., f{(Q x -2) = ©! - 1, f{(O x - 1) = 0. 
Then Vj G [[©i - 1, ©i 1]], we define /(j) =0. 

Then with the next ©2 integers [[©1 + N — s, ©1 + N — s + G>i — 1]] we construct the 
canonical ©2-cycle. We keep going on constructing for all j e [[3,/?]] the canonical 
©/-cycle. 

In consequence, we have found a function / belonging to the class [©1 , . . . , (Op}. 
Assume there exists a function g £ J^jv belonging to the class of / such that g < f. 
Let I = {i e [[0,N - 1]] such that f(i) ^ 0}. As g < f, there exists i e I such that 
■?0o) < /( z o)- There exists also jo such that io e © ;o . If /(*o) = io, then © 7o = 1, 
g(io) < io and then g(io) ^ ©/ - Then the global orbit pattern of g doesn't contain 
anymore 1 as cycle. The global orbit pattern of g is different from the global orbit 
pattern of /. If /(io) = i'o + 1» then g(io) < io- Either g(io) = io and then the global orbit 
pattern of g is changed, or g(io) < io and we are in the same situation as previously. 
Thus, in any case, the smallest function belonging to the class [©1 , . . . , C0 p ] is the one 
constructed in the first part of the proof. 

The proof of the theorem gives an algorithm of construction of the threshold function 
associated to a given gop. 

The threshold function associated to the gop [2^, l,3]n is explained in Table |9j Its rank 
is n = 25,938,474,637. 



TABLE 9. Al; 


gorithm for the threshold function construction for the ; 


gop [25,1, 3] n . 




First step 


Second step 


Third step 


Fourth step 


Fifth step 


Construction 


of Construction of 


Construction of 


Construction of 


Filling 


the re- 


the first canoni- the last canonical 


the canonical 1- 


the canonical 2- 


maining 


; images 


cal 2-cycle 


3-cycle 


cycle 


cycle 




with 




-> 1 


-► 1 


-► 1 


-» 


1 


- 


-> 1 


1 


1 -> 


1 -> 


1 -t 





1 - 


-> 


2 


2 -> 


2 -> 


2 -t 




2 - 


-> 


3 


3 -> 


3 -> 


3 -» 




3 - 


-> 


4 -> 


4 -> 


4 -> 


4 -> 




4 - 


-> 


5 -> 


5 -> 


5 -> 


5 -) 


6 


5 - 


-> 6 


6 


6 -> 


6 -> 


6 -> 


5 


6 - 


-> 5 


7 


7 -> 


7^7 


7 -> 


7 


7 - 


-> 7 


8 


8^9 


8^9 


8 -> 


9 


8 - 


-> 9 


9 -> 


9 ^ 10 


9 ^ 10 


9 -» 


10 


9 - 


-4 10 


10 -> 


10 ^ 8 


10 ^ 8 


10 — 5 


• 8 


10 - 


-> 8 



Theorem 3 There are exactly 2 W — 1 different global orbit patterns in J^y. 



That is 



(10) 



For example, for N = 4, ^(^4) = 2 4 - 1 = 15. 

Proof Let p an integer between 1 and N. Consider the set L{p,N) of p-tuples 
(a h ...,a p ) e (W) p such that ct\ + . . . + a p < N. 

We write L(N) = {L(p,N),p = 1...N}. L(N) and Sf(^)v) nave the same elements. 
Then 

p=N p=N / \t \ 

m*N) = 1 mp,n) = e „ =2^-1. 
p=i p=i \ f / 



Ordering the global orbit patterns 
We define an order relation on 5f (j^jv)- 

Proposition 1 Let A and 5 be two global orbit patterns of £f(^v). 
We define the relation -< on the set (^#) by 

A^fl iffTr(A)<Tr(B) 

Then the set (^(^n), -<) is totally ordered. 

Proof As the order -< refers to the natural order of N, the proof is obvious. 

Let r > 1, p > 1 be two integers. Let [o)i, . . . , ft) p ] and [&){,..., ft)/] be two global orbit 
patterns of £f (j^)v)- F° r example, if p <r, in order to compare them, we admit that we 
can fill [coi , . . . , (O p \ with r — p zeros and write [©1 , . . . , ft)^] = [ft)i , . . . , ft) p , 0, . . . , 0] . 

Proposition 2 Let r > 1, p > 1 be two integers such that p < r. Let A = [©1, . . . , Op] 
and B = [co[ , . . . , co' r ] be two global orbit patterns. 

. If r = p = 1 and ft>i < ft){ then A -<B. 
• If r > 2 then 

* If ©1 < ft){ then A -< 5. 

* If (0i = co[ then there exists K e [[2;r]] such that c% 7^ (0^ and Vi < i£, 
6) ; = ft)/. 

• If |A| < |5|,thenA -<B. 

• If |A| = then if co K < co' K then A -< B. 

For example, for N = 5, the global orbit patterns are in increasing order : [1] -< [1^] -< 
[I5] -< [1,2] -< [1 3 ] -< [1 5 ,2] -< [1,2,1] -< [1,3] -< [1 5 ] -< [I3, 2] -< [1 5 ,2,1] -< [1 5 ,3] -< 



[1,2, 1 5 ] -< [1,2 5 ] -< [1,3,1] -< [1,4] -< [2] -< [2,1] -< [2,1 5 ] -< [2 5 ] -< [2, l^] -< [2,1,2] -< 
[2 5 , 1] -< [2,3] -< [3] -< [3, 1] -< [3, 1 5 ] -< [3,2] -< [4] -< [4, 1] -< [5]. 

Algorithm for ordering the global orbit patterns : a pseudo-decimal order 

The Table [101 gives a method for ordering the gop : indeed, we consider each gop as if 
each one represents a decimal number : we begin to order them in considering the first 
order (0\. Considering two gops A = [a>i, cd p ] and A' = [a>[,...,a)' r ], if (0\ <(£>[, then 
A -< A'. For example, [2, 1, 2] -< [4, 1]. If co\ = a>[ and |A| - Oh < \A'\ - co[, then A -< A'. 
For example to compare the gop [1,2] and the gop [1^], we say that the first order (0\ 
stands for the unit digit - which is C0i = 1 here, then the decimal digits are respectively 
0.2 and 0.111. We calculate for each of them the modulus-C0i : we find |[1,2]| — 1 =2 
and |[1 3 ]| -1 = 3, thus [1,2] -< [1 ? ]. Finally, if (»i = Co[ and \A\ - 0} x = \A'\ - G)[, 
then also we use the order of the decimal part. For example, [l<r] -< [1,2, lj because 
1.1111 < 1.211. Applying this process, we have the sequence of the ordered gop for 
N = 4- given in the previous paragraph. 



TABLE 10. Ordered gop for N = 5 with modulus and modulus-0)i 



Gop 


Modulus 


Modulus- C0\ 


Gop 


Modulus 


Modulus-O)i 


[1] 


1 





[2] 


2 





N 


2 


1 


[2,1] 


3 


1 


N 


3 


2 


[2,1 5 ] 


4 


2 


[1.2] 


3 


2 


[2 5 ] 


4 


2 


[14] 


4 


3 


[2.I3] 


5 


3 




4 


3 


[2,1,2] 


5 


3 


[1,2,1] 


4 


3 


[2 5 ,1] 


5 


3 


[1,3] 


4 


3 


[2,3] 


5 


3 


[I5] 


5 


4 








[13,2] 


5 


4 


[3] 


3 





[1 5 > 2 ,1] 


5 


4 


[3,1] 


4 


1 


[12,3] 


5 


4 


[3,1 5] 


5 


2 


[1,2,15] 


5 


4 


[3,2] 


5 


2 


[1,2 5 ] 


5 


4 








[1,3,1] 


5 


4 


[4] 


4 





[1,4] 


5 


4 


[4,1] 


5 


1 








[5] 


5 


1 



For example, for N = 5,we construct one branch of a tree with (0\ = 1 (see Fig. [3]) : 
each vertex is an ordered orbit, the modulus of the gop is written on the last edge. 
However, the sequence of ordered gop differs from the natural downward lecture of the 
tree and has to be done following the algorithm. 




5[1, 1,1,1,1] 



FIGURE 3. Branch of the tree for the construction of the gop with C0\ = 1 on ^(#5) 

CARDINAL OF CLASSES 

In this section we emphasize some closed formulas giving the cardinal of classes of 
gop. Recalling first the already known formula for the class [l^jyv f° r which we give a 
detailed proof, we consider the case were the class possesses exactly one &-cycle, the 



case with only two cycles belonging to the class and finally the main general formula of 
any cycles with any length. We give rigorous proof of all. The general formula is very 
interesting in the sense that even using computer network it is impossible to check every 
function of J^jy when N is larger than 100. 



Discrete maps with 1 -cycle only 

The theorem 4 gives the number of discrete maps of J^y which have only fixed points 
and no cycles of length greater than one. This formula is explicit in [2D and |9(]. A 
complete proof is given here in detail. 

Theorem 4 Let k be an integer between 1 and N. The number of functions whose 
global orbit pattern is (i.e. belonging to the class [1^]a0 is ( J N N ~ k . 



That is 

h mi . , _ / 

N-k 



f N \ 

Proof Let k be a non-zero integer. Let / be a function of ^m- There are I ^ 1 

possibilities to choose k fixed points. There remain N — k points. Let p be an integer 
between 1 and N — k. We assume that p points are directly connected to the k fixed 
points. For each of them, there are k manners to choose one fixed point. There are kP 
ways to connect directly p points to k fixed points. There remains N — k — p points that 
we must connect to the p points. There are tftlpW-zt functions. Finally, the number of 

N\ N - k 



functions with k fixed points is I ) ^ kP %\\p\N-k- We now prove recursively on N 

for every 0<k<N that $[l~\ N = * J N N ~ k . We have = 1 . The formula is 

true. 

We suppose that Wk < N ft [T^ = ( NN ~ k - 

Let X be a set with N+l elements. We look for the functions of ^n+i which have k 
fixed points. Thanks to the previous reasoning, we have 

N+l-k 



We use the recursion assumption. 

N+l-k 



m^=( N _ N k +i) {N+1)N ~ Mq - e - d - 

Discrete maps with fc-cycle 

We look now for the number of functions with exactly one fc-cycle. 

Theorem 5 Let k be an integer between 1 and N. The number of functions whose 
global orbit pattern is [k]^ is Ufl^jiv x (k—l)\. 

mN = ^Nx(k-l)\. (12) 

Proof There are (Tj ways of choosing k elements among N. Then, there are 
(k— 1)! choices for the image of those k elements in order to constitute a fc-cycle 
by /. We must now count the number of ways of connecting directly or not the re- 
maining N — k elements to the fc-cycle. We established already this number which is 

equal to £ kP ^ N - k . Finally, we have ft = (*- 1)!($ ^ N - k . That is, 

P =l P=i 

tt[% = tt&x(*- 1 )U-e.d. 



Discrete maps with only two cycles 

We give the number of functions with only two cycles. 

Theorem 6 Let N > 2. Let p and q be two non-zero integers such that p + q < N. 
Then, 

HM* = HF^U = (A ,_ (p+ , ))! , • d3) 



Proof We consider a function / which belongs to the class [1^^]at- We search the 
number of functions constructed from / whose gop is [p, q]^. From the p fixed points of 



/, we construct a />-cycle. Thus, there are ( p ^ ) q _ l ) ways to choose p — l integers among 
the p + q — 1 fixed points. Counting the first given fixed point of /, we have p points 
which allow to construct (p — l)! functions with a p-cycle. Then there remain q points 
which give (q — 1) ! different functions with a g-cycle. Finally, the number of functions 

whose gop is [p,q] N is : (^j 1 ) (p - 1) \{q - 1) ! that is the formula , 



Remark We notice that for all non-zero integer such that k < N — 1, 



General case : discrete maps with cycles of any length 

We introduce now the main theorem of the section which gives the number of gop of 
discrete maps thanks to a closed formula. 

Given a global orbit pattern or, the next theorem gives a formula which gives the 
number of functions which belong to a. 

Theorem 7 Let p > 2 be an integer. Let [coi, . . . , (0 p ]n be a gop of ^(^v). Then, 

. . . , a> p ] N = m^njN^ x ( ^_| +^" x + ^ x ^+ . ..+«*,) (14) 

Ufa, ... , = ^ ^ — p (15) 

(#-(©! + . .. + ©,))! n(L ©/) 

Jfc=2 y=* 



Proof We consider a function / which belongs to [1 ;p-q: ffl ]w- We search the 
number of functions constructed from / whose gop is [o?i, . . . , (Q p ]n- From the ft>i fixed 
points of /, we construct a 0)\ -cycle. Thus, there are (° >1+ C3 ^\'~ ) ways to choose 
©i — 1 integers among the + . . . + (o p — 1 fixed points. Counting the first given fixed 
point of /, we have (0\ points which allow to construct ((D\ — 1)! functions with a 
(0\ -cycle. Then, the first fixed point of / which has not be chosen for the 0)\ -cycle, will 
belong to the ©2-cycle. Thus, there are ( C ° 2+ 0) ^ p ~ ) ways to choose ft>2 — 1 integers 
among the 0D2 + ■ ■ ■ + co p — 1 fixed points. So we have 0D2 points which allow to construct 
(fi>2 — 1) ! functions with a tt)2-cycle. We keep going on that way until there remain (0 p 
fixed points which allow to construct {(O p — 1)! functions with a (O p -cycle. Finally, we 
have constructed : 

r V-r 1 ) ^ - 1 ) r^tr 1 ) - 1) 'x.-.x r^v)K-i-i)!(%-i)! 

functions. We simplify and obtain the formula. 



Corollary 1 Let p be a non-zero integer. Let \(0\ , . . . , (o p ]n be a gop of £f (J^iv) • We 
suppose that there exists 7 such that COj > 2. Let /z be an integer between 1 and (Oj — 1. 
Then 



tt[©i, . . . , (o h . . . , (o p ] N = tt[<oi, ...,a>j-h,h, o)j+u . . . , to^ x + + • • • + <o P ). 

(16) 



Proo/ tt[»i, • • • , a>j ~ Kh, (o j+ i,. . . , (O p ] N x(h+ + ... + a> p ) = ttfl^^^JV 

(0)l+... + (Dp-l)!(/i+a)y + l+...+(Dp) 

We simplify and we exactly obtain 



tt[®i, • • • , (Qj - h,h, (o j+h . . . , (o p } N x(h+ a j+ i + . . . + (Op) = |t[fi)i, . . . , (Oj, . . . , Op]^. 
Examples : 



tt[2 5 ,l,3] n = 11,180,400. 

tt[5,2, 10,8, 15,2,3] 50 = 29, 775, 702, 147, 667, 389, 21 8, 762, 343, 520, 975, 006, 
348, 329, 578, 044, 480, 000, 000, 000, 000, 000. 

#[5,2,10,8, 15, 2, 3] 50 = 2.98 x 10 63 among the 8.88 x 10 84 functions of ^ 50 . 

FUNCTIONS WITH LOCAL PROPERTIES 
Locally rigid functions 

Obviously it is not possible to transpose to the functions on finite sets the notions of 
continuity and derivability which play a dramatic role in mathematical analysis since 
several centuries. In fact the class ^bCO of the continuous functions on the real interval 
/ is a very small subset of the set 7 R of all the functions on /. Hence by analogy to this 
fact and trying to mimic some others properties of continuous functions, we introduce 
some subsets of particular functions of J£)v, which have local properties such as locally 
bounded range in a sense we precise further. Limiting the range of the function in a 
neighbourhood of any point of the interval induces a kind of "rigidity" of the function, 
hence we call these functions locally rigid functions. In these subsets, the gop are 
found to be fully efficient in order to describe very precisely the dynamics of the orbits. 
We first consider the very simple subset Jz?.^i. N of functions for which the difference 
between f(p) and f(p + 1) is drastically bounded. In next subsection we consider more 



sophisticated subsets. 
We consider the set : 

SfXijt={fe& N such that Vi?,0<p <N-2,\f(p)-f(p + l)\ < 1}. 



Orbits ofjf&ijf 

Theorem 8 If / G Jzf^ijv then / has only periodic orbits of order 1 or 2. 

Proof We suppose that / G J^Mi^ has a 3-cycle. We denote (a;f(a);f 2 (a)) 
taking a the smallest value of the 3-cycle. If a < f(a) < f 2 (a) then there exist 
two non-zero integers e and e' such that f(a) = a + e and f 2 (a) = f(a) + e'. Thus, 
f 2 (a) - e' < f(a) < f 2 (a) + e'. That is f(a) <a< f(a) + 2e'. And finally we have the 
relation a + e <a which is impossible. 

If a < f 2 (a) < f(a) then there exist two non-zero integers e and e' such that 
f 2 (a) = a + e and f(a) = f 2 (a) + e' . Thus, f(a) - e < f 3 (a) < f(a) + e. That is 
f{a) — e < a < f(a) + e. But f(a) — e = a + e'. And finally we have the relation 
a + e' <a which is impossible. 

We can prove in the same way that the function / can't have either 3-cycle or greater 
order cycle than 3. 



Numerical results and conjectures 

We have done numerical studies of the ^(JfMi^) for N = I to 16, using the brute 
force of a desktop computer (i.e. checking every function belonging to these sets). 

The Tables [TTj [T2[ [T3| [141 H3 and [T6l show the sequences for Jzf« u to Jz^ u6 . 

In theses Tables we display in the first column all the gop of ^(Jzf^uv) for every 
value of N. For a given N, there are two columns; the left one displays the cardinal of 
every existing class of gop (- stands for non existing gop). Instead the second shows 
more regularity, displaying on the row of the gop [2r] the sum of the cardinals of all the 
classes of the gop of the form [2,2, , 1 2] which exist. 

V v ' 

k+l orders 

Then we are able to formulate some statements which have not yet been proved. 
Statement 1 



TABLE 11. Numbering the locally rigid functions for / G Jz?,«u, / G 
V.* 



g.o.p. N=l 


N=l 


N=2 


N=2 


N=3 


N=3 


N=4 


N=4 


Total number 


1 




4 




17 




68 


[1] 1 


+ 


2 


+ 


7 


+ 


26 


+ 


N 


+ 


1 


+ 


4 


+ 


14 


+ 


[%] 


+ 




+ 


1 


+ 


4 


+ 


[i 3 ] 


+ 




+ 




+ 


1 


+ 


[2] 


+ 


1 


1 


4 


4 


18 


18 


[2,1] 


+ 




+ 


1 


1 


3 


4 


[1,2] 


+ 




+ 




+ 


1 


+ 




+ 




+ 




+ 


1 


1 



TABLE 12. Numbering the locally rigid functions for / G 

■ <■ , -- f G ^i,6, / G _Sf^!, 7 . 



g.o.p. 


N=5 


N=5 


N=6 


N=6 


N=7 


N=7 


Total number 




259 




950 




387 


[1] 


95 


+ 


340 


+ 


1,193 


+ 


[1 5 ] 


50 


+ 


174 


+ 


600 


+ 


[I3] 


16 


+ 


58 


+ 


204 


+ 


[1?] 


4 


+ 


16 


+ 


60 


+ 




1 


+ 


4 


+ 


16 


+ 






+ 


1 


+ 


4 


+ 


[17] 




+ 




+ 


1 


+ 


[2] 


70 


70 


264 


264 


952 


952 


[2,1] 


12 


18 


45 


70 


166 


264 


[1,2] 


6 


+ 


25 


+ 


98 


+ 


[2 5 1 


4 


4 


18 


18 


70 


70 


1] 


1 


1 


4 


4 


17 


18 


[1,2 5 ] 




+ 




+ 


1 


+ 


[2,1,2] 




+ 




+ 




+ 


[2 5 ] 




+ 


1 


1 


4 


4 


[2 5 ,1] 




+ 




+ 


1 


1 



Statement 2 

_ N 

tt[2fc]jSf« liW = W^tfL&Sliw for k<~. 



Statement 3 

tt[2=W ljJV = tf[V for 2* < < 3* - 1. 



TABLE 13. Numbering the locally rigid functions for / e Jzf.*i,8, / G 



g.o.p. 


N=8 


N=8 


N=9 


N=9 


N=10 


N=10 


Total number 




11 Q 1 A 
1 l,o!4 




Af\ <f\1 




1 j,oy4o 


[1] 


/i 1 1 a 


+ 


1 a r\r\i 
14,U01 


+ 


4/, 064 


+ 


n 1 


mo 


+ 


/; Q<o 


+ 


OO Qf\A 


+ 


n l 


or\r\ 
/00 


+ 


1/;/; 
Z,j66 


+ 




+ 


n l 
[I4J 


1 /i 
z!4 


+ 


o/i 
/4Z 


+ 




+ 




60 


+ 


216 


+ 


754 


+ 


ri l 


16 


+ 


60 


+ 


ZlO 


+ 


ri l 


4 


+ 


16 


+ 


60 


+ 


r 1 i 


1 


+ 


4 


+ 


16 


+ 


r 1 i 

[15] 


~ 


+ 


1 


+ 


4 


+ 


r 1 1 

[2] 




+ 




+ 


1 


+ 


3,356 


3,356 


11,580 


11,580 


39,364 


39,364 


1 1 

[2,1] 


590 


952 


2,062 


3,356 


7,072 


11,580 


h oi 
[l,2j 


j6z 


+ 


i,zy4 


+ 


4,j0o 


+ 


[25] 


264 


264 


952 


952 


3,356 


3,356 


ro 1 1 


^0 
62 


/O 


000 

222 


0£ A 

264 


/ /o 


nco 
952 


ri n 1 

[!> 2 2] 


6 


+ 


28 


+ 


113 


+ 


[2,1,2] 


2 


+ 


14 


+ 


69 


+ 




1 


1 O 

to 


/u 


/O 


O/; A 

264 


0£ A 

264 


ro~ 1 1 


A 

H 


A 


lo 


lo 


ov 


/u 


1 2~ 

L ' 3J 




+ 




+ 


1 


+ 


[2,1,22] 




+ 




+ 




+ 


[25,1,2] 




+ 




+ 




+ 


[2 3 ] 


1 


l 


4 


4 


18 


18 


[24,1] 




+ 


1 


1 


4 


4 


[1,2 3 ] 




+ 




+ 




+ 


[2,1,23] 




+ 




+ 




+ 


[25.I.25] 




+ 




+ 




+ 


[25,1,2] 




+ 




+ 




+ 


I 2 *! 




+ 




+ 


1 


1 



Statement 4 



= £|f[2,2, ,^,...,2]^ ljV for2fc+l</Y 

i=1 ( th 



V 

Jt+1 orders 



ifc+i 

£ H[2pj, l,2~]^ ljV for 2k+ 1 < iV (20) 

i=l 



TABLE 14. Numbering the locally rigid functions for / E Jz?.s? U i, / £ Jf,« U2 , / € 

-2^1,13. 



g.o.p. 


N=l 1 


N=l 1 


N=12 


N=12 


N=13 


N=13 


lotal number 




45 /, /95 




1 ci c no£ 
1,515,926 




4,9 /9,/ / / 


[i] 


1 jo,oz9 


+ 


CU OA A 

J 10,o44 


+ 


1 AO'S AT2 


+ 




7S 9Q9 


+ 


OAf, 7AO 
Z^tO, / oz 


+ 


OAT 


+ 




9n OQP, 


_i_ 

1 


OJ,JJU 


T 


978 S80 
/ o,jou 


T 


[15] 




+ 


07 QfkL 


+ 


Q1 ARR 
y l,^foo 


+ 


[1-1 


Z, / 00 


+ 




+ 


OR 71R 
Zo, / jo 


+ 


n 1 

6 


9<\A 


+ 




+ 


R TiM 
0, / jU 


+ 


II7J 


9 1 A 
Z10 


+ 


7SA 

/ JO 


+ 


z,jyz 


+ 


n 1 
[%] 


01) 


+ 


Z16 


+ 


/56 


+ 


n 1 
[I9] 




+ 


aa 
60 


+ 


Z16 


+ 


n 1 

r 1 1 


4 


+ 


16 


+ 


60 


+ 


1 


+ 


4 


+ 


16 


+ 


r 1 i 


~ 


+ 


1 


+ 


4 


+ 


r 1 i 

173] 

[2] 




+ 




+ 


1 


+ 


132,104 


132,104 


438,846 


438,846 


1,445,258 


-1 A A c r\ c 

1,445,258 


[2,1] 


23,941 


39,364 


80,108 


132,104 


265,548 


438,846 


M ol 

[1,2] 


15,423 


+ 


51,996 


+ 


173,298 


+ 




11,580 


11,580 


39,364 


39,364 


1 A/I 

132,104 


132,104 


P5, 1] 


2,634 


3,356 


8,883 


1 1,580 


OA /CCA 

29,659 


OA O £. A 

39,364 




429 


+ 


1 C C C 

1,555 


+ 


5,478 


+ 


[2,1,2J 


293 


+ 


1 1 /II 

1,142 


+ 


4,227 


+ 


ro 1 

P3] 


OCT 

9jZ 


/lei 

952 


^ ^ ca 
3,356 


^ ^ c^ 
3,356 


1 1,580 


11 C O A 

1 1,580 


ro~ 1 1 


9^ 

ZjJ 


ZO^f 




CKO 

yjz 




^ ^A 
J, J JO 


ri 1 i 

[1,23] 


7 


+ 


35 


+ 


1 CO 

152 


+ 


[2,1,25] 


z 


+ 


16 


+ 


06 


+ 


[2j, 1,2] 




+ 


Z 


+ 


OA 

zu 


+ 


ro 1 


7A 


10 


OA/1 
Z04 


OA/I 
Z04 


QCO 

yjz 


QCO 

yjz 


[27,1] 
L 4 ' J 


18 


18 


70 


70 


263 


264 


[1,2tJ 




+ 




+ 


1 


+ 


[2,1,23] 




+ 




+ 




+ 


[22.1.%] 




+ 




+ 




+ 


[23,1,2] 




+ 




+ 




+ 


P , 


4 


4 


18 


18 


70 


70 


[25,1] 


1 


1 


4 


4 


18 


18 


[L25] 




+ 




+ 




+ 






+ 


1 


1 


4 


4 


[ 2 6.1] 




+ 




+ 


1 


1 



Statement 5 



1 if k = 1 

' w+l 1 (^)(Jfc+l)x3* for3<fc<^±i 



TABLE 15. Numbering the locally rigid functions for / e JSf^w, 
/ S?*i,i S . 



g.o.p. 


N=14 


N=14 


N=15 


N=15 


'1', . 1 , , 1 

lotal number 




1 O /I AO A 

16,246,924 




52,694,573 


h l 


ceil 010 


+ 


17 0/11 0/I7 
1 /,S41,z4/ 


+ 


n l 


9 A.09 9^8 
Z,0U9,Z9S 


+ 


9Q1 9A.0 


+ 


lh\ 




+ 


9 QOzL ^09 

Z,7U4,9VZ 


+ 


[i~i 


zy / , / zo 


+ 


yoz,sss 


+ 


fi~l 


94,440 




307 848 




[1-1 


29,050 


+ 


95,676 


+ 


[1-1 


S, / T-U 


1 

T 




1 

T 


[1-1 


9 ^Q9 


+ 


8 7zL8 
S, / 4S 


+ 


[I9J 


7SA 
/ JO 


+ 


9 ^Q9 


+ 




zlo 


+ 


/JO 


+ 




aa 
60 


+ 


1 A 

z!6 


+ 




to 


+ 


60 


+ 




4 


+ 


16 


+ 


ri 1 


1 


+ 


4 


+ 






+ 


1 


+ 


[2] 


4,725,220 


4,725,220 


15,352,392 


15,352,392 


[2,1] 


873,149 


1 A A C i"\ C O 

1,445,258 


2,851,350 


+ 


[I, 2 ] 


572,109 


+ 


1 DTI OIA 

1,873,870 


+ 




/I 1 O O AC 

43s, 846 


A 1 O O A C 

438,846 


1 A A C ICO 

1,445,258 


1 A A C C O 

1,445,258 




9s, 135 


100 1 a a 

132,104 


oOO oi A 

322,310 


438,846 


[1,25] 


18,873 


+ 


/CO A /; ~7 

63,967 


+ 


[2,1,2] 


15,096 


+ 


52,569 


+ 


ro~i 
[2 3 J 


9Q 9A/I 


9Q 9A/I 
jy, 904 


1 90 1 M/1 
19Z, 1U4 


1 90 1 M/1 
19Z,1U4 


ro 1 1 

i 2 3, 1J 


1U,40U 


1 1 son 


9/1 0/1 ^ 
j4,S4j 


9Q 9A/1 


[1,23] 


AC 

605 


+ 


000 

2,282 


+ 


[2,1,25] 


isy 


+ 


1 COA 


+ 


[ 2 2, 1,2] 


1 OA 

1/6 


+ 


A/t 1 

641 


+ 


[24J 


9 9CA 
9,9J0 


9 9CA 
9,9j0 


1 1 con. 

1 1,JSU 


1 1 con 

1 1,JSU 


V>~ 1 1 


QzL9 
:?4Z 


CK9 

V9Z 


^ 9Q9 


T. 9^A 

9,990 


n o 1 


Q 




+ 


/I /I 
44 


+ 


[2,1,23] 


2 


+ 


18 


+ 






+ 


2 


+ 


[23,1,2] 




+ 




+ 


P5] 


264 


264 


952 


952 


[25. 1] 


70 


70 


264 


264 






+ 




+ 


Pb] 


18 


18 


70 


70 


[2^1] 


4 


4 


18 


18 


[2,] 


1 


1 


4 


4 


[2 7 ,1] 




+ 


1 


1 



Remark We call = V[)^E^\\^^\n- F° r ^ > 2, then is the sequence A120926 
On-line Encyclopedia of integer Sequences : it is the number of sequences where is 
isolated in ternary words of length written with {0, 1,2}. 



These statements show that first the set Jzf^ijv is an interesting set to be considered 
for dynamical systems and secondly the gop are fruitful in this study. However the set 

&3hfl={fz&N such that Vp, < p < N - 2, \f(p) - f(p + 1) | < 2} 

is too much large to give comparable results. Then we introduce more sophisticated sets 
we call sets with locally bounded range which more or less correspond to an analogue 
of the discrete convolution product of the local variation of / with a compact support 
function 0%. 



Orbits and patterns of locally rigid function sets 

Consider now the set : 

„=Ue*« such that Vp, < p < iV - r - 1 , E ? a r \f{p) - f(p + r)\ < q} f]{f G 

r=t 

such that Vp,t < p < N - 1, £ OC r \f(p) ~ f(p ~ r )\ < ^} f° r the vector c£ = 
(ai,a 2 , e for<7 e N. 



TABLE 17. 


Numerical study of the set 


^ ^d t ,q.N 


for N = 10, t = 5, 


OJ] = 20, a2 = 


9,a 3 = 5, 


(X4 = 2 and as = 


1, for<7 = 


: 20, ...,142 


q maximal period 


modulus gop 


number 


functions number 


20 


1 


1 


1 


10 


26 


2 


2 


3 


82 


44 


2 


3 


6 


21,764 


49 


3 


3 


7 


48,112 


50 


3 


3 


7 


53,210 


56 


3 


4 


9 


208,692 


59 


4 


4 


15 


330,800 


63 


4 


5 


19 


626,890 


66 


4 


10 


37 


952,228 


67 


4 


10 


46 


1,064,316 


72 


5 


10 


50 


1,630,018 


74 


6 


10 


60 


1,816,826 


76 


6 


10 


61 


2,152,450 


77 


6 


10 


88 


2,416,368 


78 


6 


10 


91 


2,762,434 


79 


6 


10 


97 


3,188,080 


80 


6 


10 


99 


3,735,666 


84 


6 


10 


100 


5,876,324 


85 


6 


10 


103 


6,473,288 


87 


6 


10 


105 


7,851,728 


88 


7 


10 


121 


8,644,178 


89 


8 


10 


129 


9,521,920 


91 


8 


10 


136 


11,414,556 


92 


8 


10 


165 


12,454,440 


94 


8 


10 


175 


14,756,058 








Following next page 



TABLE 17. (Next) 



q 


IllaAlIIlal JJCllUU 


m r\ni l In c 
II1UUU1US 


gUJJ ilUIIlUCl 


tunftinnc nn m nAi* 
lUIlVUUIlo 11U111UC1 


95 


8 


10 


177 


16,077,780 


96 


8 


10 


184 


17,208,654 


97 


8 


10 


185 


18,369,854 


98 


8 


10 


188 


19,585,746 


100 


8 


10 


192 


22,083,852 


101 


8 


10 


199 


23,584,452 


102 


8 


10 


204 


25,513,892 


103 


8 


10 


244 


27,912,772 


104 


8 


10 


304 


30,560,238 


105 


9 


10 


333 


33,516,466 


106 


9 


10 


380 


36,682,960 


107 


9 


10 


424 


40,004,280 


108 


10 


10 


491 


43,685,352 


109 


10 


10 


517 


47,655,856 


110 


10 


10 


529 


51,785,410 


111 


10 


10 


562 


55,907,120 


112 


10 


10 


583 


60,341,276 


113 


10 


10 


612 


64,930,790 


114 


10 


10 


647 


69,766,178 


115 


10 


10 


706 


74,989,752 


116 


10 


10 


747 


80,087,120 


117 


10 


10 


791 


85,570,272 


118 


10 


10 


820 


91,206,218 


119 


10 


10 


836 


97,040,288 


120 


10 


10 


852 


103,121,916 


121 


10 


10 


872 


109,650,464 


122 


10 


10 


896 


116,345,296 


123 


10 


10 


919 


123,241,156 


124 


10 


10 


924 


130,360,938 


125 


10 


10 


928 


137,636,628 


126 


10 


10 


930 


145,536,068 


127 


10 


10 


932 


154,370,862 


128 


10 


10 


938 


164,145,928 


129 


1 A 
10 


1 A 

10 


960 


1 /4,94z,026 


130 


10 


10 


986 


186,438,038 


131 


10 


10 


1,006 


198,594,118 


132 


10 


10 


1,013 


211,550,402 


133 


10 


10 


1,015 


225,324,700 


134 


10 


10 


1,021 


239,976,118 


135 


10 


10 


1,022 


255,106,866 


137 


10 


10 


1,023 


286,726,234 


142 


10 


10 


1,023 


374,355,356 



The functions belonging to these sets show a kind of "rigidity": the less is q, the more 
"rigid" is the function, this "rigidity" being modulated by the vector ot t . Furthermore, 
the maximal length of a periodic orbit increases with q, and so the number of gop 
W{££s%—± ) and the maximal modulus of the gop. 

OCf ,q,N 



Remark Using this generalized notation, one has : Jzf.^i 5 „ = J£@l\^ n and 

As an example, we explore numerically the case : N = 10, t = 5, OL\ = 20, (X2 = 9, 
OC3 = 5, OC4 = 2 and OC5 = 1, for q = 20, . . . , 142. The results are displayed in Table [T7I 
In this Table "modulus" means the maximal modulus of the gop belonging to this set 
for the corresponding value of q in the row, "gop number" stands for ^(JfM— » ) and 

CCf ,q,N 

"functions number" for \iJ>?3?^ . One can point out that for the particular function 

CCf ,q,N 

ot t of the example; it is possible to find 10 intervals I\,l2,---,ho C N such that if 
q E I r then there is no periodic orbit whose period is strictly greater than r, (e.g., 
h = [[74,87]]). Furthermore it is possible to split these intervals into subintervals I rs in 

which f > ) is constant when q thumbs L f . This is not the case for tLSf^ 1 — > 

* V CC t ,q.Nj * r ' S * a t .q,N 



CONCLUSION 



A discrete dynamical system associated to a function on finite ordered set X can only 
exhibit periodic orbits. However the number of the periods and the length of each are not 
easily predictable. We formalise such a gop as the ordered set of periods when the initial 
value thumbs X in the increasing order. We can predict by means of closed formulas, 
the number of gop of the set of all the function from X to itself. We also explore, using 
the brute force of computers, some subsets of locally rigid functions on X, for which 
interesting patterns of periodic orbits are found. Further study is needed to understand 
the behaviour of dynamical systems associated to functions belonging to these sets. 
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TABLE 16. Numbering the locally rigid 
functions for / e _£?« 116 . 



g.o.p. 


N=16 


N=16 


Total number 




1 '~if\ r\f^ , —ir\r\ 

170,028,792 


[11 
I 1 ] 


S7 4.77 S49 


_i_ 
i 


1-1 


26 912 398 




1~ 


9 114 088 


+ 


1~ 


3 097 650 


+ 




77O, /0+ 


+ 


n l 


1 /I c/: 


+ 




96,096 


+ 




29,158 


+ 


[19] 


8,748 


+ 


r 1 1 

[%] 


2,592 


+ 


756 


+ 


[%] 


216 


+ 


[%] 


60 


+ 




16 


+ 


[%] 


4 


+ 




1 


+ 


[2] 


49,610,818 


49,610,818 


[2,1] 


9,255,822 


15,352,392 


[1,2] 


6,096,570 


+ 


H 


4,725,220 


4,725,220 


[25, 1] 


1,051,686 


1 A A c r\ c 

1,445,258 


r 1 ^ i 

[1,2 5 ] 


213,975 


+ 


[2,1,2] 


179,597 


+ 


\H 


438,846 


438,846 


[25, 1] 


114,798 


132,104 


r 1 ^ i 

[1,25] 


8,284 


+ 


[2,1,25] 


6,146 


+ 


[25,1,2] 


2,876 


+ 


n 1 






[2~ 11 


11,246 


11,580 


[1,2 3 ] 


204 


+ 


[2,1,23] 


106 


+ 


[22,1,25] 


22 


+ 


[23,1,2] 


2 


+ 


H 


3,356 


3,356 


[2 5 ,1] 


951 


952 


[L25] 


1 


+ 


[2gl 


264 


264 


[26,1] 


70 


70 


[27] 


18 


18 


[27.I] 


4 


4 


[2«] 


1 


1 



